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Abstract 

Quantum-classical escape-rate transition has been studied for two general 
forms of magnetic anisotropy in ferrimagnetic or antiferromagnetic particles. 
It is found that the range of the first-order transition is greatly reduced as 
the system becomes ferrimagnetic and there is no first-order transition in al- 
most compensated antiferromagnetic particles. These features can be tested 
experimentally in nanomagnets like molecular magnets. 
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Understanding the mechanism of transition between two states in a bistable system is 
of utmost importance of the fast-growing area of macroscopic quantum phenomena (MQP). 
|jlj] The rf SQUID and the closely related current-biased Josephson junction |2j have been 
typical examples of MQP. In recent years, as another good candidates for MQP, small 
magnetic particles have been investigated via quantum tunneling of magnetization. [[J In 
general, transition between two states can occur due to two mechanisms. At sufficiently 
high temperature the transition rate T obeys the Arrhenius laws, T ~ exp(— ?7//cbT), with 
U being the height of the energy barrier. At a temperature low enough to ignore the thermal 
fluctuation quantum tunneling is relevant with T ~ exp(—U/hu) where u is the oscillation 
frequency in the well. Thus, the crossover temperature T is expected to exist between the 
thermal activation and quantum tunneling. Whether the two regimes smoothly join or not 
has been intensive subject over the past few years. j§-|§ 

The criterion of first- or second-order transitions for the escape rate has been suggested 
by Chudnovsky, [Q who showed that if the oscillation period r(E) of the instanton is mono- 
tonic (nonmonotonic) function where E is the energy of the instanton, dT{T)/dT becomes 
continuous (discontinuous), i.e., the second- (the first-) order transition around the crossover 
temperature. Since then, Gorokhov and Blatter || carried out the nonlinear perturbation 
near the top of the barrier, obtained a criterion for the first-order transition in tunneling 
problem based on Ref. ||, and applied it to the phase transition of vortex creep pinned 
in the columnar defect. Later, theoretical investigations for the transition in spin systems 
have been performed by several groups. |6|-|9| Until now theoretical studies have been fo- 
cused on the ferromagnetic particle. However, considering that most ferromagnetic systems 
are actually ferrimagnetic particles, the strong exchange interaction should be taken into 
consideration. The exchange interactions are expected to suppress the first-order transition, 
and thereby, if it is large, there may be no first-order transition. Thus, it is clearly desirable 
to assess its importance for the phase transition before contemplating real experiment. In- 
deed, the exchange interaction in the Hamiltonian does not permit a simple mapping onto 
the particle problem ||§ and the periodic instantons in the reduced one-dimension. |7||| 
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None of such previous theoretical methods are applicable to the present spin system with 
the exchange interaction. In this paper we will develop new approach to treat the phase 
transition in ferrimagnetic or antiferromagnetic particles based on the spin-coherent-state 
path integral, and present complete analytic results of the phase boundary between first- and 
second-order transitions for two general forms of the magnetic anisotropy energy. Also, we 
will discuss effects of the fourth order anisotropy in the Hamiltonian on the phase transition 



in molecular clusters. 10 



Let us consider a small ferrimagnetic or antiferromagetic particle with two magnetic sub- 
lattices whose magnetizations, mi and ni2, are coupled by the strong exchange interaction 
m! •m 2 /x_L, where x± is the perpendicular succeptibility. In the case of a non-compensation 



of sublattices with m(= mi — m 2 > 0), the Euclidean action can be taken to be JTT 



S E (9,<f>) = V Jdr [I ™ 1 ± m2 4> - i™<j> cos 9 + C 1E j , (1) 

where V is a volume of the particle, 7 the gyromagnetic ratio, and 9, 4> spherical coordinates 
of sublattice magnetization mi which determines the direction of the Neel vector. Here, 
C\e is the Euclidean Lagrangian density which is composed of the exchange interaction, the 
anisotropy energy and the energy given by an external magnetic field. The first term in 
the action (Jl|) known as the topological term JT2J generates the phase factor related with 



quantum phase interference. In this work we do not include this term because it is not the 
purpose of this paper to discuss this issue. 

In the simplest easy-axis anisotropy along the z-axis and a transverse magnetic field 
along the x-axis the Euclidean Lagrangian density is expressed as 



C 1E = £±[(0 + iryH sin 0) 2 + U sin 9 + i^H cos 9 cos 0) 2 ] 
27^ 

+K\\ sin 2 9 — mH sin 9 cos (/>, (2) 

where x± = XiX m il m \)- Then, the corresponding classical trajectory is determined by the 
equations 

— in u 4> sin 9 + x u [9 — -<fi 2 sin 29 + 2ibcf) cos sin 2 9 



-i& 2 sin2#cos 2 0] - 2cos#(sin# - /icos0) = 0, (3) 
in u 9 + x u [(p sin 9 + 20$ cos 9 — 2*66* cos sin # 

+ ife 2 sin20sin#] - 2/isin0 = 0. (4) 

where we have introduced m/(if||7) = n u , X±/(-^||7 2 ) — ^u, jH = b, H/H c = h, and 
H c = 2K\\/m. Note that H c is the critical magnetic field at which the barrier vanishes. In 
this case we decompose 9 {(f)) into the position of the top of the barrier 9 (0) and a fluctuation 
term 59 (50), i.e., 9 = 9 + 59 (0 = + 5(f)) for the behavior of the weakly time-dependent 
solutions, where 9 = n/2 and = 0. Our goal is to solve Eqs. (||) and fl4|) for 59 (r) and 
50(r) and find the correction to the oscillation period away from the thermal saddle point. 
|J Denoting 5f2(r) = (59 (r), 50(r)), we have 5f2(r + /?/i) = 5fl(r) at finite temperature and 
write it as Fourier series 5fl(r) = J2™=-oo exp(zu; n r) where a; n = 2Ttn/f3Ti. Since simple 
analysis shows that 59 is real and 50 imaginary in this model, to lowerst order we write 
them in the form 59 ~ a9 u i cos(a;r) and 5<f> ~ ia<f> u i sin(u;r) in the vicinity of the thermal 
saddle-point solution. Here we have parametrized the solutions of the equation of motion by 
the amplitude a of the oscillations, which quantifies the difference between the thermal and 
the time-dependent solutions near the top of the barrier. Substituting them into Eqs. ([|) 
and (II), we obtain the relation 



Ki Xu\u l u± - V) - 2 + 2h -n u u u ± 



Qui n u u u ± x u (uj 2 ± — b 2 ) + 2h' 

where n u = n u — 2bx u , and the oscillation frequency 

u2 n 2 u -2(l-2h)x u 



(5) 



"" ± = ° 2^ 



u 



[hi - 2(1 - 2h)x u Y + 16A(1 - h)xl - m\xl 
- ^ ' <6) 

In order to see the behavior of the oscillation period via the frequency we need to investi- 
gate Eqs. (^|) and (f|) by writing 59 ~ a9 u \ cos(co>r)+56 l 2 and 50 ~ iacf> u i sin(a;r)-|-i502, where 
592 and 502 are 0(a 2 ). Neglecting terms of order higher than a 2 , we find u = u> u+ , 56*2 = 0, 
and 502 = 0. This implies that there is no shift in the oscillation frequency. Thus, to third 



order in perturbation theory we use 56 ~ a6 u \ cos (cur) + 56% and 5 (ft ~ ia4> u i sin (cur) + i5(p 3l 
where 56 3 and 5cf> 3 are 0(a 3 ). Substituting them into Eqs. @ and ([|), and neglecting terms 
of order higher than a 3 , we obtain for the change of the oscillation frequency 

niyl^-ul + )^~u^)=a 2d fg u (h,y u ). (7) 

Here y u = x u /n 2 l (= x±K\\/m 2 ) indicates the relative magnitude of the noncompensation. 



For large noncompensation (y u <C 1, i.e., m ^> Jx±K\\) and for small noncompensation 



(y u ^> 1, i.e., m <C yJx±K\\), the system becomes ferromagnetic and nearly compensated 
antiferromagnetic, respectively. Without any approximation, g u {h,y) in Eq. ([?]) is so com- 
plicated that it is not illuminating to present its specific form. Noting that the first-order 
region is expected to shrink by increasing y u , it is slightly influenced by the magnetic field 
terms in the exchange interaction of Eq. (0). Thus, neglecting those terms in the exchange 
interaction, we can obtain the analytic form of g u (h,y) approximaterly given by 

1 

? 



g u {h, y)~-[2-2y + 12hy - 3y 2 + 8hy 2 - 2y 3 



-(2 + 2y + Ahy - y 2 )^l - Ay + Shy + Ay 2 }. (8) 

As shown by Chudnovsky, [Q if the oscillation period r is not a monotonic function of a 
where a is a function of E in the absence of dissipation, the system exhibits a first-order 
transition. Thus, the period r(= 2-n/uj) in Eq. (0) should be less than r u+ {= 27r/u; M+ ), 
i.e., u) > u u+ for the first-order transiton. It implies that g u (h,y u ) > in Eq. (f7|) for the 
first-order transition, and g u (h,y u ) = determines the phase boundary between the first- 
and the second-order transition. Using Eq. (||]), we approximaterly get the phase boundary 

{Vu) - 6 12/(i/ u )V» + I2y u ' W 



where 



f( y ) = -18y 2 - 350y 3 - 372y 4 + 66y 5 - y e 



+6 V3y 3 (l + 3y + y 2 ) 2 (16 - 13y + 32% 2 - 6y 3 ). (10) 



Also, in order to compare Eq. (§) with the phase boundary from the exact calculation, 
we have employed the numerical method for the exact phase boundary, whose results are 
illustrated in Fig. [I]. The approximate result is found to be slightly different from the 
exact one obtained by the numerical methods. The scaled magnetic field, h for the phase 
boundary decreases linearly for y u <C 1, as h ~ 1/4 — y u and parabolically for y u < 1/2, 
as h ~ 0.397(y u — 1/2) 2 . Thus, it is evident that the region for the first order transition is 
greatly reduced as the system becomes ferrimagnetic and there is no first-order transition 
in almost compensated antiferromagnetic particles. 

We now turn to the biaxial symmetry with the Euclidean Lagrangian density 



C 1E = ^[9 2 + (4>sm 



2 7 2 

+Ki cos 2 9 + K 2 sin 2 9 sin 2 0, (11) 



where K\ > K 2 > are the anisotropy constants. This model describes a hard axis z and 
an easy axis x, and the corresponding classical trajectory satisfies 

— irib(j) sin 9 + Xb 0) — -4> 2 sin 29^j 

+ sin 2^(1 - fcsin 2 0) = 0, (12) 
in b 9 + Xb (j) sin 9 + 2(j)9 cos 9^j 

— 2fcsin#sin</>cos0 = 0. (13) 

where we have introduced m/(Kij) = rib, X±/(-^i7 2 ) = x b, and K 2 /Ki = k. In the high 
temperature regime the solution of Eqs. (|1^) and (|13| ) is 9 = = tt/2. In order to find the 
order of the transition, let us expand Eqs. ( |T2"{ ) and ( |T3"D into a series around this solution, as 
9 = it/2 + 89 and = 7r/2 + <50. Simple analysis for Eqs. (|Tl~|) -(|r3"l) shows that 89 is imaginary 
and 8cj) real. Thus, to lowest order in perturbation theory, we write 89 ~ ia9b\ sin(a;r) and 
8(j) ~ a(pbi cos(ujt). Substituting them into Eqs. (0) and (p~3|) while neglecting terms of 
order higher than a, we have 

_ x b Ub± + 2(1 - k) _ n b uj b ± 
9 b i n b uJb± x b ujl ± -2k 



and the oscillation frequency 

n\ + 2(1 - 2k)x h 



2x1 



v ni + 4(l-2k)n 2 b x b + Ax 2 b 
± 2x1 ' ( } 

Next, let us write 88 ~ %a8 b \ sin (cur) + i88 2 , and 50 ~ a<f> b \ cos{ur) + 502, where 5#2 and 

502 are of the order of a 2 . Inserting them into Eqs. ( |12"D and (0), we arrive at uo = uo b+ and 

882 = 502 = 0. In order to find the change of the oscillation period, we proceed to the third 

order of perturbation theory by writing 88 ~ %aB\,\ sin(cG>r)+z5#3, and 50 ~ a<f) b i cos (cur) +503. 

Substituting them again into Eqs. ( p~2|) and ( p~3|) and retaining only terms up to 0(a 3 ), we 

have 

<vl^ 2 ~ <)(^ 2 - = a 2 ^g b (k,y b ), (16) 

where yb = Xb/n b (= x±Ki/m 2 ) and 

2 



^(fc, 2/) = -Jl + 3y - Aky - Ay 

y 
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- (1 + y - 2y 2 )^l + 4y- Sky + Ay 2 }. (17) 

Here we note that the parameter y b has the same physical meaning as y u in uniaxial symmetry 
by replacing K\\ by K\. Continuing in the present case as in the earlier model we have 
g b {k,y b ) > for the first-order transition and the phase boundary given by 

k= l -{l-y b ){l + 2y b ) 2 . (18) 

As follows from Fig. 0, the criterion f?| that in case of the biaxial symmetry the first-order 
transition occurs at k > 1/2 is not valid in ferrimagnetic particles. Fig. ^| shows that the 
ratio of two anisotropy constants, k increases linearly for y b <C 1, as k ~ (1 + 3y b )/2 and 
has a maximum at y b = 1/2. Therefore, the region for the first-order transition is greatly 
suppressed with increasing y b and vanishes beyond yb = 1/2. 

To illustrate the above results with concrete examples we discuss two molecular clusters, 
Mn 12 and Fe 8 . Actually, both samples are ferrimagnetic, and thereby y u and y b should be 



taken into account in uniaxial and biaxial symmetry, respectively. Taking the measured value 
of the anisotropy parameter, e.g., k = 0.728 for Fe§, it is seen from Fig. |^ that the system 
may exhibit the first-order transition for yj, < 0.157 and the second-order transition beyond 
that region. This point should be checked in experiment. Also, since neutron scattering data 
[TIJ shows that an Fe§ (Mn 12 ) should be described by the biaxial (uniaxial) spin Hamiltonian 
with the fourth order anisotropy term, it is meaningful to discuss how much it makes an effect 
on the phase transition. Our analysis shows that the first-order region is slightly reduced 
by adding the fourth order terms in the Hamiltonian. Taking Fes for instance, the range 
of yb for the first-order transition is changed to < yj, ^ 0.153, as shown in Fig. |3|. Thus, 
considering the fact that most ferromagnetic systems are ferrimagnetic, it is very important 
to obtain the information on the magnitude of the quantities y u and y^ for observing the 
first-order transition in real experiments. 

In conclusion, we have presented the phase diagrams for first- and second-order tran- 
sitions in ferrimagnetic or antiferromagnetic particles. It is found that the range of the 
first-order transition is greatly suppresed by the exchange enhancement. In fact, this can 
be qualitatively understood from the consideration that the exchange interaction produces 
some effective magnetic field and thereby plays an important role in reducing the range of 
the first-order transition. In this respect, the first-order transition can not be observed in 
almost compensated antiferromagnetic particles, which is consistent with our analytic re- 
sults. This general features is expected to be observable in nanomagnets including molecular 
clusters. 

Discussions with D. A. Gorokhov, W. Wernsdorfer and B. Barbara are acknowledged. 
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FIGURES 

FIG. 1. Phase diagram h(= H/H c ) vs y u (= X^K\\/m 2 ) obtained by the exact calculation for 
the uniaxial symmetry in a transverse magnetic field. Note that h and y u contain information of the 
magnetic field and the relative noncompensation, respectively, and thereby the system tends to be 
ferrimagnetic or antiferromagnetic with increasing y u . See the text for details. Inset: Comparison 
of the exact result (a) with the approximate one (b). 

FIG. 2. Phase diagram k(= K2/K1) vs yb(= X^Ki/m 2 ) for the transition in biaxial symmetry. 
Note that k contains information of the relative strength of the biaxial two axes, and the first-order 
transition is expected for k < 1 in which the system starts to be biaxial from uniaxial. Also, the 
parameter yb has the same physical meaning as y u in uniaxial symmetry. See the text for details. 

FIG. 3. Change of the range of the first-order transition for the biaxial spin Hamiltonian in 
Fes (a) with and (b) without the fourth order anisotropy term in the Hamiltonian. Note that (a) 
< yb ^ 0.153 and (b) < yb < 0.157 for the first-order transition. 
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